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Abstract: Torque ripple occurs in Permanent Magnet Synchronous Motors (PMSMs) due 
to the non- sinusoidal flux density distribution around the air-gap and variable magnetic 
reluctance of the air-gap due to the stator slots distribution. These torque ripples change 
periodically with rotor position and are apparent as speed variations, which degrade the 
PMSM drive performance, particularly at low speeds, because of low inertial filtering. 
In this paper, a new self-tuning algorithm is developed for determining the Fourier Series 
Controller coefficients with the aim of reducing the torque ripple in a PMSM, thus 
allowing for a smoother operation. This algorithm adjusts the controller parameters based 
on the component's harmonic distortion in time domain of the compensation signal. 
Experimental evaluation is performed on a DSP-controlled PMSM evaluation platform. 
Test results obtained validate the effectiveness of the proposed self-tuning algorithm, with 
the Fourier series expansion scheme, in reducing the torque ripple. 
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1. Introduction 

The Permanent Magnet Synchronous Motor (PMSM) market is growing more rapidly when compared 
to traditional competitors because of lower cost, as well as higher efficiency and reliability. For the sake 
of energy savings and environmental performance, PMSMs also feature one of the highest torque to loss 
ratios. These motors are widely used in fast dynamic positioning systems and machine-tool 
components [1,2]. The main disadvantage of PMSMs is the non-uniformity in the developed torque, 
known as "torque ripple" [3]. Torque ripple generates speed oscillations which cause system 
performance deterioration and in machine-tool applications, it can leave visible patterns on high precision 
machined surfaces [4]. Under the assumption of pure sinusoidal back electromagnetic force (EMF), the 
conventional Field Oriented Control (FOC) applies constant current references in the synchronous 
reference frame to produce a constant torque. However, depending on the magnet shape and how well 
the windings are manufactured, the back-EMF has in practice very different waveforms, which range 
from almost sinusoidal to trapezoidal. Torque ripple occurs in PMSMs due to non- sinusoidal flux 
density distribution around the air-gap and the variable magnetic reluctance of the air-gap due to stator 
slots distribution. These torque ripples change periodically with rotor position and are apparent as 
speed variations, which degrades the PMSM drive performance, particularly at low speeds because of 
low inertial filtering [5]. 

In order to improve the performance of PMSMs and increase its market share, the suppression of 
the pulsating torque has received much attention in recent years [6-17]. These torque ripple reduction 
techniques can be divided into two groups: one focusing on the improvement of motor design and the 
other emphasizing the use of active control of stator current excitation. From the motor design point of 
view, skewing the stator lamination stacks or rotor magnets, arranging proper winding distributions 
and incorporating other motor design features reduces cogging torque partially, but does not 
completely eliminate it. Moreover, special machine design processes add addition complexity to the 
production process, which results in higher machine cost [18]. 

The second approach, which is of our interest, concentrates on using an additional control effort to 
compensate for the periodic torque pulsations. Some methods rely on pre-programed stator current 
excitation to cancel torque harmonics. However, accurate information about the PMSM parameters is 
required, and a small error or variation in these parameters can produce higher torque ripple due to 
the open-loop control. As an alternative, closed-loop control algorithms with online estimation of 
parameters and adaptive control algorithms have been proposed to reduce torque ripple. One possible 
approach relies on a closed-loop speed regulator to attenuate indirectly torque pulsations since all 
possible sources of torque ripple are observable from rotor speed, and hence this method has potential 
for complete torque ripple minimization. Repetitive Control techniques incorporate a sinusoidal 
control component to deal with periodic torque pulsations [19-23] while Iterative Learning 
Control (ILC) is implemented in the frequency domain to reduce torque ripple, by means of Fourier 
series expansion [24-29]. Some recent papers deal with learning control algorithms for Permanent 
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Magnet Step Motors [30-32], by identifying the Fourier coefficients of any truncated approximation 
and implementing Iterative Learning techniques, providing an experimental comparison for both 
methods. Adaptive techniques have been proposed based on the spectrum of the torque perturbation 
using some theoretical developments. 

In this paper, a new self-tuning algorithm is developed for determining the Fourier Series Controller 
coefficients with the aim of reducing the torque ripple in a PMSM, thus allowing for a smoother 
operation. This algorithm adjusts the controller parameters based on the component's harmonic 
distortion in the time domain of the compensation signal. The estimated Fourier coefficients are used 
by a nonlinear controller which achieves accurate and ripple-reduced torque control. Experimental 
evaluation was performed on a DSP-controlled PMSM evaluation platform and test results obtained 
verify the effectiveness of proposed self-tuning algorithm, with the Fourier series expansion scheme, 
in reducing the torque ripple. 

This paper is organized as follows: A model of the Permanent Magnet Synchronous Motor is 
presented in Section 2. The new self-tuning Fourier Coefficient Algorithm is introduced in Section 3. 
Section 4 describes the experimental setup, and the experimental results are presented in Section 5. 
Finally, in Section 6 concluding remarks are provided. 

2. Model of PMSM 



In this section, a standard PMSM model [33] is revised and additional considerations are explained 
so torque ripple sources are clarified. For a three-phase PMSM, the flux linkages ^abcs related to the 
mutual and self-inductances Ls and currents / are given as: 

^abcs = L s i(ibcs (1) 



In matrix form: 



^as ~ ^ as as v as 



i nq + L as ft s ifo s + L ascs i cs + ^asm 



^bs Lfr sas i as + Lbsbs^-bs ^bscs^cs ^bsm 
^cs ~ Lcsasi-as ^csbs^bs ^cscs^cs ^ csm 



(2) 



The stator windings voltages u abcs depend on the winding resistance r s and flux linkages babes'- 

. d*P a bcs (3) 



^abes ~ 1**1 



s L abcs 



dt 



Rewriting this expression in matrix form: 



"Uas" 




~r s 


0 


o- 




las 
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r s 


0 




i-bs 


U CS . 




.0 


0 


r s . 




- l cs- 



dt 
d x ¥ bs 



dt 

d¥ r , 



(4) 



dt 

The stator windings are displaced by 120°, and the flux linkages ^Vasm-^bsm-^csm established by 
the permanent magnet, which are periodic functions of 9 r , are assumed to be sinusoidal with 
magnitude 
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^ asm sin(0 r ) 

^bsm - sin (Or ~ \ n ) 

^csm - sin (Or + \ n ) 

From Equation (3), assuming L s is constant: 



(5) 



u afrcs — r s^abcs L t 



(6) 



dt 



dt 



d6 r 



Defining co r = — r - , we have: 



Hence, in Cauchy form, by using L s 1 : 

dlabcs 



o) r cos(0 r ) 
a) r cos (o r — 

0) r COS (j) r + ^7r) 



(7) 



dt 



d*¥ 

_ _w -1 - r-1 771 I r-1 

— ~L>s ^s^abcs ~ L> s ~~ ^ ^s ^abcs 



(8) 



Incorporating the transient behavior of the mechanical system, where electric torque T e , load torque 
T L , viscous friction coefficient B m and inertia moment / are used: 

d 2 9 rm (9) 



T e B m a) rm T L — J - ^ 2 



da) T 



dt 



-j(7e~ B m a) rm - T L ) 



(10) 



where mechanical position and velocity are related by 



dt 



= 0) T 



To find the electromagnetic torque developed T e , where W PM is the permanent magnet energy, the 
co-energy W c is used: 

¥ m sin(0 r ) 



W c — -[i as hs ics]L s 



l bs 



as l bs 



ics\ 



^mSin (6 r - ft) 
*F m sin (e r + ft) 



+ W, 



PM 



(11) 



Therefore, we have the following formula to calculate the electromagnetic torque for the 
three-phase P -pole permanent-magnet synchronous motors: 

PdWc PT™ / f _ 9 \ / . 9 \\ (12) 



Te 2 d6 r 



■ (i as cos(d r ) + i bs cos (e r - ft) + i cs cos (e r + ft)) 



Hence: 
doj 



rm r 1 in 



(i as cos(0 r ) + i bs cos (e r - ft) + i cs cos (e r + ft)) -^yu r m --Tl 



dt 



2/ 



1 

7 



(13) 



Using the electrical angular velocity a) r and displacement "0 r ", related to mechanical angular 

2 2 

velocity and displacement as 0) rm = -0) r and 0 rm = -0 r , results in the following equation: 
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da) r 
dt 



4/ 



■ (ias cos(0 r ) + i bs cos (d r - V) + i cs cos (d r + |7r)) - ^ fiV - 7} 



2/ 



(14) 



d0 r 
dt 



= 0) r 



Regarding the implicit time reparameterization to express the time functions acceleration and the 
speed in Equation (14) as functions of the rotor position, the relation rotor position- time is guaranteed 
to be invertible for the total rotor position (not only for one revolution range), if rotor position is a 
monotonic function of time like happen for a non-cero angular speed of constant sign. 

To control the angular velocity, one regulates the currents fed or voltages applied to the stator 
windings. To maximize the electromagnetic torque developed, the motor should be fed by a balanced 
three-phase current set: 

^as(0 = V2i M cos(a> r t) =V2i M cos(6) e t) =V2J M cos(0 r ) 

*cs(0 = ^2^M cos (^r^ + ^) = V2i M cos (a> e t + ^7r) = V2i M cos + ^7r) 
Generating an electromagnetic torque: 



(15) 



V2i M (cos 2 (0 r ) + cos 2 (0 r - §7r) + cos 2 (# r + §7r)) 



2 ' ™\ ^ " \ > 3 J v 3 J J 2V2 

To produce the specified current, the balanced three-phase voltages are given as: 

u as (t) = V2u M cos(0 r ) 
u bs(t) = V2u M cos (o r — ^7r) 
u cs (t) = V2u M cos (o r + ^7rj 



(16) 



(17) 



To simplify the control of PMSM, it is a common practice to transform the equations from 
three-phase voltages u abcds to the qdO variables u qd0s for the rotor reference frame [34]. 

In matrix form, the mathematical model of the PMSM in the rotor reference frame is given as: 
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(18) 



The required currents to regulate the angular velocity of PMSM and guarantee balanced operating 
conditions are given as: 



i$ s (t) = V2/ M , i r ds (t) = 0, i£ s (t) = 0 
And assuming that inductances are negligible, the applied voltages should be: 

u r qs (t) = V2u M , u r ds (t) = 0, u r Qs (t) = 0 



(19) 



(20) 
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2.1. Additional Considerations 

In order to understand torque ripple in PMSM we have to reconsider some assumptions from the 
previous model. 

2.1.1. Non-Sinusoidal Flux Linkages 

Flux linkages are not perfectly sinusoidal so the electromotive force differs from cosine function 
and applying cosine currents to the stator windings produces torque ripples. The induced non-cosine 
electromotive forces e abcs are all assumed to be periodic functions with a peak value E p [35]: 

fastfr) 
fbs (Or ~ f 1 ) 



abcs ~ ^ — £*pto r 



fcs (?r + \ n ) 



(21) 



where the functions f as (6r)> fbs(fir)> an d f C s(fir) have the same shape as e as ,e bs , and e cs with a 
maximum magnitude of ±1. 

Since electromagnetic torque is given by: 

■ ■ ■ JL ( 22 ) 

*e ~ l e as l as e bs l bs e cs l cs\ xv 

(x) r 

Using Equation (21), we finally have: 

T e = Ep [fas^r)i a s + fbs (?r ~ \ n ) l bs + fcs (?r + \ n ) ks] ^ 

where T e is still independent of frequency, but in this case current waveforms should be calculated to 
produce constant torque. 

2.1.2. Non-Constant Inductances 

From Equation (1), electromotive force can be calculated by: 

_ <M abcs _d (24) 

^abcs ~ ^ ~ dt abcs * rnJ 

Without assuming constant L s : 

digbcs dL s d*¥ m (25) 

&abcs ~ L s , + labcs j. * i. 

at at at 

The second term from Equation (25) produces torque ripple due to inductance angular variations, 
and it is associated to differences in winding inductances of the stator. 

2.1.3. Stator Yoke Reluctance Variations 

From a macroscopic viewpoint, the torque produced in a PMSM is given by [35]: 

T -^ — -^jv 2 dR (26) 
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As mentioned before, the first term appears when motor construction causes the winding inductance 
to vary as a function of position, and third term describes the mutual torque that is used to make the 
motor shaft turn. Additionally the second term describes cogging torque that appears whenever rotor 
magnetic flux travels through the varying reluctance of stator yokes, attempting to align with the stator 
teeth or poles independent of any current. When motor shaft is rotated by hand, the pulsations felt are 
caused by cogging torque. 

3. Self-Tuning Fourier Coefficient Algorithm 

As stated in the previous section, torque ripples arise from non-sinusoidal flux density distribution 
around the air-gap and variable magnetic reluctance due to stator slots distribution. These torque ripples 
change periodically with rotor position and are apparent as speed variations particularly at low speeds. 

These periodic torque ripples T r (0 r ), with period r = 2n can be represented in the form of a 
Fourier series: 

N 

T r (Q r ) = a 0 + ^[a^ cos(/c0 r ) + b k sin(/c0 r )] ^7) 

k=l 

where a 0 , a k , and b k are unknown constant vectors. 

Since cogging torque and harmonic components of the non-sinusoidal electromotive force depend 
on the slot distribution, torque ripple is a periodic function of the position and can be considered 
anti-symmetric and modeled by the sinusoidal components. 

" (28) 
T r (Q r ) = 2^[h sin(fc0 r )] 

k=l 

Considering the inertia moment of the system /, the acceleration a r , angular velocity 0) r , in relation 
to 0 r , are given by: 

a r (0 r ) = = y[^[b k sin(k0 r )] 



(N 
^[-^cos(k0 r ) 



(29) 



(30) 



vk=l 

To compensate for the velocity ripple a control voltage u(t) should be introduced where K v is the 
voltage constant: 

Kjsr\b 



N - ' (31) 



U(0 r ) = -y( ^ |-^COS(/C0 r ) 



Introducing c k = the control voltage can be written, as function of the position 0 r (t), as: 

Ik 

N 



^(9r) = ^[^cos(/c0 r (t))] (32) 

k=l 
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For each term c /c cos(/c0 r (t)) the angular position ripple Q r i PP (t) can be approximated by a sinusoidal 
function, and because of speed variations its temporal representation of cos(/c0 r (t)) is distorted, thus 
we can use a measure of its distortion S to iteratively adjust its coefficients c k until distortion is 
reduced by regulating o) r (0 r ). Let r = 2% be the spatial period for a complete mechanical revolution: 

(33) 



r T 

5= J cos (kQ ripp (t) ) dt 
Jo 



c k (t) = c k (t-l) + 8S (34) 

This algorithm permits adjusting the control voltage parameters, adapting for changes in the torque 
ripple, and the parameter 5 allows for controlling the adjusting speed. 

Demonstration: 

While the shape of the cogging torque is a complex function of motor geometry and material 
properties, here it is approximated by a sinusoidal function and consequently the angular position 
ripple is also approximated by a sinusoidal function. 

Assuming that angular position is given by 0 r (t) = a) 0 t + Asin(a) 0 t + cp), were nominal angular 
position "a) 0 t" comes from the nominal constant speed of the rotor "a) 0 '\ and the angular position 
ripple "Asin(a) 0 t + cp)" corresponds to the angular lag or advance (from nominal angular position), 
produced by periodic perturbation such as cogging torque. Then the angular velocity 0) r , and angular 
position ripple Q ripp are given by: 

6ri W (0 = Asin(a) 0 t + cp) (35) 

a) r (t) = a) 0 + Aa) 0 cos((jL) 0 t + cp) (36) 

For each term c /c cos(/c0 r (t)), according to Fourier series methodology, its value is in theory 
determined by: 

- T (37) 



1 f 

Ck = j J cos (k0 ripp (t)) co r (t)dt 



were T represent the temporal period for a complete mechanical revolution. 

In practice, 0) r {t) is not necessarily a periodic function of time and synchronization of the 
controller to the torque perturbation could be difficult to achieve. 

Using Equations (35) and (36) in Equation (37): 

1 f T (38) 
c k = - I cos(k Asin((o 0 t + (p)) (a) 0 + Aco 0 cos(co 0 t + (p))dt 
T Jo 

So: 

a) 0 f T 1 f T (39) 

c k = — I cos(kAsin(a) 0 t + (p)) dt + — I Aa) 0 cos(a) 0 t + cp) cos(kAsin(a) 0 t + (p)) dt 
* Jo * J 0 

Recognizing that the second term is equal to zero, because can be rewritten as — J Q cos(u) du: 

(40) 



co 0 f 

c k = — I cos(kAsin(co 0 t + cp)) dt 
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Finally: 



Ck = yJ cos(ke ripp (t))dt 



(41) 



when the rotor speed is not constant or periodic of unknown period it is preferable to replace T by x to 
allow for a fix parameter. Furthermore instead of computing ^7, because the angular velocity is 

changing, the integral term can be used to iteratively adjust the coefficient c, through the gain factor 5 
as proposed in Equations (33) and (34). 

4. Experimental Implementation 

Figure 1 shows the overall torque ripple minimization scheme. During the transient state, the 
Fourier Series Controller is not activated and the Field Oriented Control (FOC) sets the motor in a 
stable operation. 

Figure 1. Block diagram of the Fourier series controller applied to the Field Oriented Control. 
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When steady state is reached, the Fourier Series Controller is applied and it provides the additional 
compensation so as to minimize torque ripple. Conventional PI current controllers that generate the 
control voltages in accordance with the field oriented control are used in the inner loop. The current 
controllers work with a sample time of 500 |is and gains are set as: K p = 1, Ki = 80, all variables are 
considered in per unit values and the 8 parameter is set to 0.02. 

Figure 2 shows the configuration of the experiments. A TMDSHVMTRPF development system 
with a F28035 DSP control card is connected to the EMJ-04APB22 four pole pairs Permanent Magnet 
Synchronous Motor with the following parameters: 200 V, 2.7 A, 400 W maximum power, 300 rpm 
rated speed, 4.7 Q, 0.014 H stator resistance and inductance, and 2500 PPR incremental encoder 
attached. Measurements of the variables were taken at the PWMDAC ports of the TMDSHVMTRPF 
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development system by using a FLUKE 199C floated oscilloscope-meter, setting the bandwidth to 10 
kHz for high frequency rejection. The encoder signals were coupled with a TTL Buffer (SN74LS243N 
Bus Transceiver), to avoid electric noise that degrades the angular position signal readings. 

Figure 2. Experimental setup including the TMDSHVMTRPF development and 
EMJ-04APB22 PMSM. 




The performance evaluation of the controller with the proposed self-tuning algorithm is presented in 
the following section. 

5. Experimental Results 

To verify the performance of the proposed self-tuning algorithm with the Fourier series expansion 
scheme, experiments were performed using the setup described in the preceding section. The 
experiments were conducted for speeds lower than 10% of the motor's nominal speed. The 
performance criterion used to evaluate the performance of the proposed scheme for torque ripple 
minimization is the variation of the angular speed determined from the angular position measurements 
from the encoder. 

The Standard Field Oriented Control 

Figure 3 shows the angular position of the motor versus time, the slope of the triangular waveform 
change reflecting speed variations. In Figure 4 a cosine function of the angular position, with the same 
frequency as the perturbation fundamental frequency, is plotted against time, showing a distortion 
shape caused by speed variation that can be used to adjust the Fourier coefficients. Figure 5 presents 
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the speed ripple, although it is only 1.6% of the reference speed (4.4 rpm of 273 rpm), it is responsible 
for distorting the previous two waveforms. 

Figure 3. Angular position of the motor, at 273 rpm, controlled by Field Oriented Control. 
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Figure 4. Fundamental cosine term, distorted by speed ripple, for Field Oriented Control. 
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Figure 5. Speed ripple, at 273 rpm, controlled by Field Oriented Control. 
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The proposed self-tuning algorithm with the Fourier series expansion scheme (for the first two 
terms) Figure 6 shows a reduction in angular position waveform distortion. Figure 7 verifies that the 
fundamental cosine function shape also improved. Figure 8 presents a 2 rpm peak to peak speed ripple 
which is half of the original variation as shown in Figure 5. The waveform in Figure 9 corresponds to 
the control signal u(0) in per unit values. 

Figure 6. Angular position of the motor, at 273 rpm, controlled by the proposed scheme, 
for the first two terms. 
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Figure 8. Speed ripple, at 273 rpm, controlled by the proposed scheme, for the first two terms. 
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Figure 9. Control signal u(0), of the proposed scheme, for the first two terms. 
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Figures 10-12 shows similar results for the proposed control for the first four terms. Figure 13 
presents the control signal u(0), although the wave forms in Figures 9 and 13 are substantially 
different, for this experiment the results are similar and perhaps it is because torque ripple is almost 
sinusoidal. Figures 14-17 shows that the second and fourth harmonic cosine function shape also 
improved. Finally Figure 18 shows that speed can be reduced to 2.7% of nominal speed while the 
motor starts to malfunction at 5.7% of nominal speed with only FOC control. 

Figure 10. Angular position of the motor, at 273 rpm, controlled by the proposed scheme, 
for the first four terms. 
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Figure 11. Fundamental cosine term, corrected by the proposed scheme, for the first 
four terms. 
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Figure 12. Speed ripple, at 273 rpm, controlled by the proposed scheme, for the first 
four terms. 



i i 1 1 1 i 1- 



H M M h 



275 

270 

? 
& 

s 



0 0.1 0.2 Time(s) 

Figure 13. Control signal u(0), of the proposed scheme, for the first four terms. 
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Figure 14. Second harmonic cosine term, distorted by speed ripple, for Field 
Oriented Control. 
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Figure 15. Fourth harmonic consine term, distorted by speed ripple, for Field 
Oriented Control. 
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Figure 16. Second harmonic cosine term, corrected by the proposed schem. 
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Figure 17. Fourth harmonic consine term, corrected by the proposed scheme. 




Figure 18. Angular position of the motor, at 80 rpm, controlled by the proposed scheme. 
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6. Conclusions 

In this paper has been presented an adaptive self-tuning algorithm for determining the Fourier 
coefficients of the controller with the aim of reducing the torque ripple in a PMSM. Its implementation 
is simple and represents a good alternative for minimizing torque ripple, cogging torque and non- 
sinusoidal electromotive torque variations due to its periodic nature. The proposed scheme does not 
require previous knowledge of the motor parameters. The performance of the proposed scheme has 
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been evaluated through experimentation and test results confirm 50% speed ripple reduction (from 4.4 
rpm to 2 rpm peak to peak speed ripple). Further research should be conducted to extend these results 
to applications where load torque varies as periodic function, which is not considered in this work. 
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